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SUMMARY 

Lutomirski' s  graphical  results  concerning  the  intensity  distribution  of  i 
focussed  laser  beam  in  the  presence  of  turbulence  have  been  confirmed  and 
extended  to  include  annular  apertures.  We  have  also  considered  segmented  aper¬ 
tures  and  present  some  indication  of  the  likely  improvement  in  focal  intensity 
when  the  segments  are  optimally  adjusted  to  counteract  turbulence. 


Accession  For 

Wi'  t  J  CiA&I 
DLC  IAB 
i  Unannounced 
i  JuctiT ic.  ticn 


Copyright 

© 

Controller  HMSO  London 
1979 


"3  _ 

Di st”i 

bnt  ■> 

t 

P.l'O.  a  A  * 

A/all  fti'd/oi 

_ 

A 

lal 

LIST  OF  CONTENTS 


1  INTRODUCTION  3 

2  FORMULATION  OF  THE  PROBLEM  (CIRCULAR  APERTURE)  3 

3  THE  EVALUATION  OF  THE  VARIOUS  INTEGRALS  6 

3.1  M(6,g,x)  7 

3.2  Ms(x,z)  7 

3.3  <I>(a,z)  10 

4  EXTENSION  TO  ANNULAR  APERTURES  11 

5  SEGMENTED  APERTURES  AND  ADAPTIVE  CONTROL  13 

6  DISCUSSION  OF  RESULTS  15 

7  CONCLUSIONS  17 

Appendix  A  Using  program  M21LASER  19 

Appendix  B  Using  program  M210PTSR  21 

Appendix  C  The  probability  density  function  for  the  separation  of 

two  points  placed  randomly  within  an  annulus  22 

Table  1  Showing  the  improvement  in  Strehl  ratio  predicted  by 

optimising  the  segmented  aperture  of  Fig  4b  24 

References  25 

Illustrations  Figures  1-7 

Report  documentation  page  inside  back  cover 


* 


O 

W 


M-C  7903 


3 


INTRODUCTION 


Lutomirski1  derived  a  formula  for  Che  mean  intensity  distribution  resulting 
from  a  finite  laser  beam  taking  turbulence  into  account.  He  considered  only  cir¬ 
cular  apertures  and  it  was  our  purpose  to  confirm  his  numerical  results  and  then 
to  extend  his  analysis  to  include  annular  and  sector  shapes.  The  former  was 
relatively  straightforward  but  the  latter  proved  more  awkward  in  that  a  simplify¬ 
ing  transformation  introduced  by  Lutomirski  could  not  then  be  used.  We  shall 
refer  extensively  to  his  paper  and  many  of  his  equations  are  reproduced  for  con¬ 
venience,  their  original  numbers  being  distinguished  by  the  prefix  L. 

2  FORMULATION  OF  THE  PROBLEM  (CIRCULAR  APERTURE) 

The  starting  point  for  our  work  was  equation  (L5)  for  the  mean  intensity 
distribution. 


<I>(£,z) 


r2  ‘  2E  *  (*i  ' 


X  Ms(r1,r2,P)UA(r])U*(r2)d2rId2r2 


Note  that  equation  (L5)  was  misprinted.  Fig  1  explains  most  of  the  geometric 
notation,  k  is  the  wave  number  (k  ■  2tt/A)  .  UA  describes  the  field  in  the 
aperture  and  includes  a  focussing  term 


UA(jj)  *  expj-  ir2(a  2  +  ikf_1)j 


(LI  4) 
(2) 


The  amplitude  distribution  is  thus  considered  to  be  Gaussian  with  standard 
deviation  a  .  Hence  a  uniform  distribution  is  obtained  by  setting  a  =  «  . 

Mg  is  the  spherical  wave  mutual  coherence  function  (MCF).  Turbulence 
enters  via  Mg  which  essentially  describes  the  loss  of  coherence  due  to  turbu¬ 
lence.  Mg  depends  on  the  distance  between  pairs  of  points  and  it  is  not 
necessary  to  consider  distances  in  excess  of  the  aperture  diameter,  D  . 

Lutomirski  introduces  the  transformation 


p-  "  5,  -  r2 


r  -  i(r.  +  r,) 


to  get  hi 8  equation  (L6). 


4 


< I  Kg. z)  ■  j  d^pMs(p,z)e  Cik/*)g.g  j  u^(r  +  Jg) 


x  U»(t  -  ie)eak,l'e-5d2r 


The  MCF  is  defined  by  equations  (L7)-(L10), 

OO 

/ 


Ma(p,z)  *  exp^-  ~ 
c 


1 

/  "o' 


1  - 


dKK<|>(K)  /  Jrt(Kpu)du 

0  0 


OD 

/ 


K$(K)dK 


where 


z'1  *  2irk2  j  Ki)>(K)dK 

0 


(4) 


(5) 


(6) 


2 

4> (K)  is,  to  quote  Lutomirski  and  Yura  ,  the  spectral  density  of  atmospheric 
indexr-of-refraction  fluctuations.  They  suggest  a  fairly  arbitrary  modification 
of  the  Kolmogorov  spectrum  to  extend  its  validity  to  large  and  small  wave  numbers, 
ie  instead  of 


<KK)  »  0.033C2K~'  1/3 
n 


they  propose 


<KK) 


(7) 


(8) 


where  and  Lq  are  the  inner  and  outer  scales  of  turbulence  where  respect¬ 

ively  viscous  damping  and  inhomogeneity  become  important.  Equation  (7)  is  there¬ 
fore  valid  in  the  so-called  inertial  subrange  Lq1  <  K  <<  .  Throughout  ourwork 

we  have  assumed  *  0.001  m  and  Lq  »  1  m  .  The  latter  is  especially  import¬ 
ant  and  is  often  taken  to  be  the  mean  distance  above  the  ground  which  perhaps 
.  .  .  2 

indicates  how  imprecise  a  quantity  it  is.  CR  characterises  the  strength  of  the 

turbulent  variations  of  refractive  index.  ■  0  represents  no  turbulence, 

.  2  -12  -2/3  n 

while  Cn  *  10  m  would  be  strong  turbulence.  Again,  this  is  clearly  a 

very  variable  quantity. 
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Putting  equation  (8)  in  (6)  and  using  the  fact  that  £  <  Lq  gives 


0.39k2C2L^/3 
n  0 


(9) 


We  shall  return  to  the  detailed  calculation  of  M  in  section  3.2  for  it  is 

s 

evident  from  equations  (4)  and  (5)  that  it  can  be  performed  independently  of  the 

main  integral.  Indeed  it  can  be  seen  from  equation  (5)  that  one  could  tabulate 

once  and  for  all  the  term  inside  the  square  brackets  as  a  function  of  p  , 

2 

given  Cn  . 

What  is  perhaps  not  so  obvious  at  this  moment  is  that  the  four-dimensional 
integral  in  equation  (4)  can  be  simplified.  Substituting  expression  (2)  in  (4) 
gives  for  the  inner  integral 


/  u0  exp 

-  •  r 

J 

2 

Writing  dr*  rdrdQ  (see  Fig  2)  and  taking  the  real  part  of  the  integral  gives 


4u: 


K 

CM 

(e)  r 

J  i 

j  exp  « 

0  0 

(r2  +  jp2) 


cos 


kpr(i  -  i) 


cos  0 


rdrde 


If  we  then  redefine  the  integration  variables  as  shown  in  Fig  2  and  write 
x  ■  p/D  ,  the  integral  over  the  area  of  overlap  of  the  two  circles  becomes 


cos~'x  1 

r 

I,(x)  *  D2Uq  J  dip  J  exp 

0  x/cos  p 

- 

— n  (u  +  2x  -  2ux  cos  d>> 

4a2 

x  cos 


nr  (i  •  ?)<u  cos  ♦  ■ x)] 


x) I udu 


Defining  further  dimensionless  quantities 


D/2a 


0  =  JkD' 


;H) 


(10) 


results  in  an  expression  identical  to  equation  (L15)  except  that  the  latter  is 
normalised 


6 


-1  , 

COS  X  1 


:j(<5,8,x)  *  D2Uq  j  d<j>  uexp|~<S2(u2  +  2x2  -  2ux  cos  <)>)J 

i|28x(u  cos  <|>  -  x)j< 


x/cos  $ 


x  cos  2gx(u  cos  <(>  -  x)ldu  .  (11) 


2  2 

Writing  d  p  =  pdpdt//  =  D  xdxdi|/  ,  equation  (4)  then  becomes 

2  1  71  / 

<I>(p,z)  =  ^2^^  2  j  dxD2x!j  (6,g,x)Ms(x,z)  J  cos|^  p  cos  i(/jd  t|/  , 


U2  j  2\2  ^ 

<I)(a,z)  =  f  xJ0(2ax)M(5,6,x)Ms(x,z)dx  , 


(12) 


where 


M(6,B,x)  =  Ii(6,8,x)/D2U^  and  a  =  ^  . 


(13) 


Equation  (12)  is  equivalent  to  Lutomirski's  equation  (LI 6)  and  is  in  prin¬ 
ciple  much  easier  to  evaluate  because  no  more  than  a  double  integral  has  to  be 
computed  at  any  stage.  In  practice,  some  caution  is  necessary  because  as 
a  and  |g|  increase,  both  Jq(2oix)  and  M(6,g,x)  become  highly  oscillatory. 

The  three  parameters  a,  g,  6  allow  three  of  the  major  features  to  be 
separately  studied 

(i)  a  the  effect  at  points  off-axis 

(ii)  6  the  effect  at  points  away  from  the  focal  plane 

(iii)  6  the  effect  of  a  Gaussian  rather  than  uniform  aperture  distribution. 

Lutomirski  presented  only  one  diagram  (half  of  Fig  LS)  in  which  more  than  one  of 
these  is  non-zero.  We  shall  see  in  section  4  that  the  extension  to  annular 
apertures  involves  the  introduction  of  a  fourth  parameter. 

3  THE  EVALUATION  OF  THE  VARIOUS  INTEGRALS 

3.  1  M(<5 ,  g,x) 

The  double  integral  of  equation  (11)  has  been  dealt  with  in  a  straight¬ 
forward  fashion  by  a  2-D  Simpson  procedure .  We  recognise  that  the  cosine  term 
in  the  integrand  really  necessitates  a  formula  such  as  Filon's  and  this  would 
be  a  worthwhile  modification  if  extensive  studies  were  planned  of  regions  well 
removed  from  the  focal  plane.  For  the  present  we  observe  that  the  largest 
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argument  occurring  in  the  cosine  term  is  6/2  when  x  ■  1  ,  ^  ■  0  ,  u  ■  1  . 

With  an  integration  step,  Au  0.05  we  can  reasonably  trust  our  simple-minded 
scheme  up  to  6  =*  20  .  We  have  also  adopted  A4>  «  0.05  (3°).  These  values  are 
empirical  and  were  selected  after  studying  the  accuracy  of  the  integration  pro¬ 
cedure  for  various  different  settings.  Such  an  investigation  was  possible 
because  <I)  can  be  obtained  analytically  on  the  axis  (Jq(2csx)  *  1)  in  free 
space  (Mg  =  1 ) . 

To  be  specific  (of  equation  (L20)), 


1 

I  ~  xM(6,0,x)dx 

6 


,  -6^/2  .  20// 
+  Ae  sm  0/4 


+  6 


(14) 


Note  therefore,  that  the  mean  intensity  at  the  focus  (0  =  0)  is 

,2 


<1 


(15) 


Hence  for  uniform  aperture  distributions  (6=0)  the  integral  in  equation  (12)  is 
the  Strehl  ratio  (SR)  which  is  why  equation  (12)  was  written  in  that  way. 

3.2  Ms(x,z) 

The  MCF  is  defined  by  equation  (5)  with  (8)  and  (9).  Lutomirski  quotes  a 
simple  expansion  (equation  (L 11))  which  is  valid  for  most  cases  of  practical 
interest.  Writing 


Mg(x, z)  =  exp $Dg(x,z)J 

where  Dg  is  called  the  wave  structure  function,  Lutomirski’ s  approximate 
expression  is 


(16) 


D  (x,z)  =  — 

S  Zq 


/D  \S/3  /D  \2 

■te)  -"(« 


(17) 


z  D  D 

Lutomirksi  and  Yura  plot  against  showing  how  it  tends  to  unity  when 

Dx  >  Lq  •  Our  investigations  suggest  that  equation  (17)  is  an  excellent  approxi¬ 
mation  (<1%  error)  for  Dx/Lq  4  0.3  .  Since  x,<  1  and  Lq  =  I  m  we  can 


A 


8 


therefore  employ  it  throughout  for  aperture  diameters  not  exceeding  30  cm. 

Note  that  we  are  unlikely  to  require  Dx/Lq  >  1  . 

We  now  explain  how  equation  (17)  is  obtained  and  corrected  when 
Dx  >  0.3Lq  .  From  equations  (5),  (6)  and  (16), 

CO  1 

Dg(x,z)  -  8w2k2z  f  K<f>(K)dK  f  -  J0(KDxu)ldu  .  (18) 

U  U  L 

Introducing  (8)  and  using  (9)  gives 


2  2 
~K  *0 


Dg(x,z) 


20z 

3z 


/  d0tL°>  /  ['  - 


z  D 
c  s 

4z 


5  7  ;yft)  i  r  ^  1 

6  J  ~~rr/6 dy  ]  ['  "  Jo(^Dxu/Lo)Jdu  * 


(19) 


where  y  -  (KLQ)2  . 

From  here  on  we  shall  assume  unity  for  the  exponential  factor  in  (19) 
which  is  valid  provided  that  the  contribution  to  the  integral  from  y  >  (Lq/JIq)2 
is  negligible.  The  y  domain  is  divided  into  two  regions  which  are  integrated 
separately.  The  break  point  is  yj  *  (LQ/Dx)  so  that  JQ  can  be  replaced  by 
a  power  series  for  y  <  yj ,  and  the  inner  integral  can  be  performed  analytically. 


ie 


1  -  JQ(w) 


iw2 .  a*2)2  +  a*2)3 . 

(2 ! ) 2  (3! )2 


J 


1  - 


Jq(wu) 


■] 


du 


I  4  1  6 

320  w  16128  w 


(20) 


Two  terms  are  sufficient  for  >99%  accuracy  when  w  *  1  .  The  integral  from 
y  ■  0  to  y  *  yj  can  then  be  evaluated 
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When  Dx  0. 3Lq  the  region  y  >  y j  must  be  treated  numerically  because 
the  series  expansion  is  not  applicable  when  y  >  1  .  We  have 


*  /  2  TTTTfnre  /  [’ 

G§)  0 


-  Jq(  /yDxu/LQ)Jdu 


Substituting  w  =  (y  +  1)  gives 


w. 


. 


dw 


f  (w) 


JQ(u)du 


1  - 


f  (w) 


(25) 


where  f(w)  =  (w  -  1)^Dx/Lq 

3116  Wj  =  ja0/Dx>2  +  lj 

(Hence  f(w()  =  I). 

The  integrand  decreases  from  1  at  w  =  0  to  =*  0.08  at  w  =  Wj  . 

Note  that  if  the  w  domain  is  divided  into  N  equal  intervals  for  the  purposes 
of  numerical  integration  then  the  largest  argument  of  the  Bessel  function  occurs 
at  the  smallest  w  and  is 


f(v)  ■  j»6/5[<VD!!)2  * 


,3/5 


[i  +  (Dx/Lq)2J 


(26) 


One  therefore  processes  the  integral  from  Wj  down  to  0  so  that  the  inner 
integral  can  be  evaluated  progressively.  The  value  of  zcDs/4z  is  finally 
obtained  by  adding  (21)  and  (25).  As  Dx  -*•  ®  ,  yj  -*•  0  so  Sj  vanishes  while 
Wj  -*■  1  ,  f(w)  -+ 90  giving  ^2  =  1  as  expected. 


3-3  <1  )  (ct,z) 

All  numerical  integrations  are  carried  out  using  Simpson's  rule.  Those 
described  in  the  previous  sections  use  a  constant  interval  which  is  written  into 
the  computer  program.  Acually,  each  range  is  divided  into  an  even  number  of 
equal  steps  each  no  wider  than  the  prescribed  value. 


M-C  7903 


M-C  7903 


1 1 


The  integral  in  equation  (12)  is  tackled  adaptively  using  QAOZA  of  the 

Harwell  Subroutine  Library.  This  is  because  of  the  oscillatory  behaviour  of 

Jn(2ax)  and  M(<$,g,x)  and  the  exponential  nature  of  M  (x,z)  resulting 
u  s 

in  an  awkward  integrand.  We  have  included  a  facility  which  provides  a  useful 

check.  D  is  computed  at  x  =  0.25  using  equation  (17)  and  if  it  exceeds  a 
s 

chosen  value  (ELIM),  the  regions  0-0.25  and  0.25-1  are  integrated  separately  and 

summed.  The  program  is  anyway  on  the  alert  for  M  -*  0  and  does  not  integrate 

s 

beyond  x'  where  D^x', z)  =  40  .  One  can  therefore  run  the  program  with 
ELIM  =  0  and  ELIM  =  40  and  compare  the  results.  The  former  would  normally 
be  expected  to  be  the  more  accurate. 

4  EXTENSION  TO  ANNULAR  APERTURES 


Annular  apertures  present  no  great  difficulty,  the  procedure  being  essen¬ 
tially  the  same  except  that  the  integration  for  M(6,  B.x)  is  over  the  area  of 
overlap  of  two  annuli.  Defining  £  to  be  the  ratio  of  inner  and  outer  diameters, 
four  cases  can  be  distinguished  (see  Fig  3). 


(a) 

(b) 

(c) 

(d) 


xijd+c) 

id  +  £)  >  X  >  4(1  -  £) 

x  | ( 1  -  £) 
x  <  5 


only  the  outer  circles  overlap 

inner  circle  overlaps  outer 

whole  of  inner  circle  is  within  both  outers 

inner  circles  overlap. 


It  can  be  seen  from  Fig  3  that  the  range(s)  of  u- integration  depends  on  the  par¬ 
ticular  choice  of  (j)  ,  x  and  £  .  For  example,  in  case  (d)  there  is  no  need 
to  carry  out  u-integration  at  all  at  small  $  while  in  case  (c)  at  small  4> 
there  would  be  two  distinct  regions  to  be  treated.  Consideration  of  the  various 
configurations  leads  to  the  following  rules 

(i)  No  integration  necessary  when  cos  $  >  max 

(ii)  When  sin  <J>  £/2x  ,  the  full  u-range  is  valid. 

(iii)  When  sin  <f>  <  £/2x  ,  solve 


/x  1  +  4x2  +  £^ 
\£  ’  4x  )  ' 


2 

v 


4vx  cos  0  +  4x 


2 


0 


for  Vj  and  v ^  (vj  K  and  integrate  over 

the  range  Vj  v<  u  <  Vj  should  it  overlap. 


u  from  x/cos  <p  to  1  except  in 


12 


The  formula  equivalent  to  (14)  is 


f 


—  xM(6,6,£,x)dx 


A 


4  (e-£252/2  -  a'62'2)2  4  4a-<'*52>S2/2  sin2^(,  .  ?2)J 


(64  +  82)  (27) 


which  is  again  valuable  for  verifying  numerical  procedures.  Note  that  when 

o  2 

5  =  6  =  0,  the  value  of  the  integral  reduces  to  (1  -  £*■)  . 

An  interesting  feature  of  the  double  integral  for  11(6, 6, x)  (or 
M(6,6,£»x))  is  that  when  6=6=0  it  reduces  to  an  expression  which  is 

proportional  to  the  area  of  overlap  of  two  identical  circles  (annuli)  and  hence 
related  to  the  probability  density  function,  P(x)  for  the  separation  of  two 
points  placed  randomly  in  a  circle  (annulus)  eg  consider  equation  (11)  with 

6  =  6  =  0, 

-1 

cos  X  l 

D2M(0,0,x)  =  J  d$  J  udu 

0  x/cos  $ 


=  ID 


and 


/  -1  /  2\ 

ICOS  X  -  XM  -  X  I 

16  (  -1  2/  FN 

=  —  I  X  COS  X  -  x  /  1  -  x  I 

*  v  / 


(28) 


xM(0,0,x) 


(29) 


Hence  for  6*6*0,  the  integral  in  equation  (12)  for  points  in  the  focal 
plane  can  be  written 


! 


j  P(x) Jq (2ax)Mg (x, f )dx 


(30) 


Gilbey  has  derived  the  same  result  (for  a  ■  0  ,  in  which  case  the  integral  is 
the  Strehl  ratio)  from  first  principles  without  following  Lutomirski's  line  of 
development. 
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Note  that  it  is  possible  to  integrate  (30)  analytically  in  free  space. 


This  is  the  well-known  Airy  pattern  for  the  focal  plane  intensity. 


(31) 


5  SEGMENTED  APERTURES  AND  ADAPTIVE  CONTROL 

The  Strehl  ratio  denotes  the  reduction  of  intensity  due  to  turbulence. 

It  is  clear  from  equations  (16)  and  (17)  that  for  given  z ,  Mg  will  be  larger 

when  D  is  small.  Thus,  to  maximise  the  Strehl  ratio  the  smallest  possible 
aperture  should  be  used  to  reduce  the  maximum  separation  of  pairs  of  points. 

However,  the  reduced  aperture  will  give  rise  to  a  reduction  of  diffraction- 

limited  intensity  at  focus  in  free  space.  This  effect  can  be  overcome  by  syn¬ 
thesising  a  large  aperture  from  a  number  of  smaller  apertures  or  'segments' 

(see  Fig  4a).  The  wave  launched  from  each  segment  must  be  independently  con¬ 
trolled  in  phase  so  that  the  contributions  from  the  various  segments  are  in 
phase  when  they  reach  the  focus.  In  other  words  the  objective  would  then  be  to 
adjust  the  phase  of  each  segment  to  enhance  the  intensity,  particularly  at  the 
focus.  In  this  approach,  it  must  be  remembered  that  the  Strehl  ratio  obtainable 
cannot  exceed  that  due  to  a  single  element,  ie  the  degradation  of  intensity  on 
account  of  turbulence  will  be  at  best  as  if  the  total  power  was  being  channelled 
through  a  single  segment.  If  the  aperture  contains  segments  of  different  shapes 
and  sizes  (see  Fig  4b)  then  the  maximum  attainable  Strehl  ratio  will  be  some 
intermediate  value.  Therefore  we  need  only  consider  one  segment  of  each  type 
and  to  combine  the  results  appropriately  in  order  to  predict  the  optimum  focal 
intensity.  Note  that  it  is  not  possible  to  calculate  the  corresponding  intensity 
at  other  points  because  the  relative  phases  of  the  aperture  segments  are  not 
actually  determined. 

Since  circular  symmetry  is  now  lost,  transformation  (3)  is  no  longer  help¬ 
ful  and  we  must  return  to  equation  (1)  (as  indeed  was  preferable  for  the  deriva¬ 
tion  of  equation  (14)  or  (27)).  Considering  only  the  focus,  the  mean  intensity 
is,  using  (2)  , 
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The  integration  limits  can  be  reduced  because  M  only  involves  the  separation 

s 

of  pairs  of  points  (see  Fig  4c)  and  not  the  direction  in  which  that  separation 
is  measured.  This  is  analogous  to  the  previous  treatment  involving  overlapping 
circles  where  only  one  quarter  of  the  region  of  overlap  needed  to  be  considered. 
Hence 
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2  2 

Introducing  x^  =  (2rj/D)  ;  x„,  *  (2r2/D)  and,  as  before,  C  *  Rj/R2  > 

6  =  D/2a  ,  gives 


<I>(f)  =  4-1 
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III 
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-<Xj+x2)5  /2 


Ms(v,f)d62dx2d01dx,  (34) 
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2  2 

where  v  =  $D  x(  +  x2  -  2/xjX_2  cos  (@2  -  0  j )  | 


It  is  then  easily  verified  that  in  free  space  (M  =  1)  , 


(I>(f) 


(35) 


so  by  setting  0  ■  2ir  the  intensity  due  to  an  annular  aperture  is  regained. 


A  computer  program  using  a  four-dimensional  Simpson  scheme  has  been  written 
to  evaluate  (34)  for  an  aperture  of  the  type  depicted  in  Fig  46  and  some  results 
are  discussed  in  the  next  section.  We  now  explain  how  the  Strehl  ratio  for  the 
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complete,  optimally  adjusted  aperture  is  determined  from  the  values  for  the 
individual  sectors. 


When  the  aperture  is  a  solid  circle  or  annulus  one  may  still  think  of  it  as 
being  made  up  of  sectors  which,  in  the  presence  of  turbulence,  are  non-optimally 
adjusted  so  that  at  the  focus,  when  the  complex  disturbances  are  summed,  the 
amplitude  of  the  sum  is  less  than  the  sum  of  the  amplitudes.  Clearly,  the  latter 
represents  the  best  possible  performance. 


Consider  an  aperture  consisting  of  K  annuli  such  that  the  outermost 

circle  is  of  diameter  D-  =  D  and  the  innermost  is  D.  and  define 

0  k 


Ei 


■  VDi-l 


=  SR 


for  any  sector  in  the  ith  annulus. 
The  optimised  Strehl  ratio  is  given  by 


(36) 


where  I^  is  the  focal  intensity  in  free  space  due  to  the  ith  annulus.  Thus 
for  the  configuration  of  Fig  4b  and  a  uniform  aperture  distribution 
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If  the  central  circle  had  also  been  illuminated,  £3  =  0 

and  -  /r^(l  -  +  '/n2^(I  "  ^2)  +  ^3^ 


opt 


6  DISCUSSION  OF  RESULTS 

Two  FORTRAN  computer  programs  have  been  written  to  evaluate  the  integrals 
in  expressions  (12)  and  (34).  The  first  has  been  used  to  verify  the  graphical 
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results  of  Lutomirski^  (Figs  L4-L8),  excellent  agreement  being  obtained.  How¬ 
ever,  it  is  to  be  expected  that  relative  accuracy  will  be  lost  as  the  Strehl 
ratio  decreases  or  as  a  and/or  8  increase  (te  as  the  integrand  in  (12) 
and/or  (11)  becomes  more  oscillatory).  This  is  not  particularly  disturbing 
since  one  is  mainly  interested  in  the  intensity  at  the  focus  and  its  immediate 
vicinity  while  we  are  confident  that,  actually  at  the  focus,  discrepancies  are 
negligible  for  n  0.001  .  It  is  unlikely  that  cases  with  n  <  0.001  have  any 
practical  significance  so  that  our  neglect  of  the  exponential  term  in  4>(k) 
(equations  (8)  and  (9))  can  be  justified  as  follows. 

Ignoring  the  inner  scale  of  turbulence  implies  that  Mg(Jl.Q,z)  =  1  and 
hence  that  separations  less  than  do  not  contribute  significantly  to  (I>  . 

Clearly,  the  region  0  4x4  £q/D  will  only  be  important  if  the  integrand  of 
(12)  decreases  rapidly  with  x  or  if  D  «  .  We  do  not  envisage  even  the 

smallest  element  of  a  segmented  aperture  having  an  incircle  less  than  lOt.^  dia¬ 
meter  while  the  first  eventuality  implies  a  very  small  Strehl  ratio.  It  can  be 
expressed  as  a  limitation  on  range  from  equation  (17) 


6.3A2/C2  (38) 

n 

-12  -2/3 

is  in  units  of  10  m 

Some  results  obtained  with  the  main  program  (see  Appendix  A  for  details  of 
its  implementation)  for  an  annular  aperture  are  presented  in  Figs  5  and  6  and 
are  compared  with  those  for  a  circular  aperture  of  the  same  outer  diameter  (see 
also  Figs  L5b  and  L6b).  Fig  7  may  be  compared  with  Fig  L8a. 


z  ^0.39k2C2A3/3j 


where  lr 


0.001  m  ,  A  is  in  pm  and  C 


The  second  program  has  been  used  to  predict  the  maximum  Strehl  ratio 

attainable  with  the  aperture  shown  in  Fig  4b,  the  results  appearing  in  Table  1. 

Appreciable  enhancement  of  the  Strehl  ratio  is  clearly  theoretically  possible 

2  -12  -2/3 

although  in  this  example  the  factor  of  15  achieved  at  *  10  m  would 

still  be  considered  insufficient  and  one  would  have  to  divide  the  same  size 

.  4 

aperture  into  many  more  segments  in  order  to  reach  90pt  ^0.5  .  Gilbey  has 
confirmed  these  calculations  using  equation  (30)  by  approximating  the  sectors 
by  rectangles,  for  which  P(x)  is  available  in  an  analytical  form.  Appendix  B 
explains  how  to  use  our  program  while  Appendix  C  contains  analytic  expressions 
for  P(x)  for  an  annulus. 


Neither  program  consumes  more  than  a  few  seconds  of  central  processor  time 
on  the  ICL  1906S  (and  run  especially  rapidly  when  D/Lq  <  0.3)  and  they  are  also 
more  accurate  than  Monte-Carlo  multi-dimensional  quadrative  procedures. 
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7  CONCLUSIONS 

Two  FORTRAN  programs  have  been  written 

(i)  to  calculate  the  mean  intensity  distribution  due  to  a  focussed  laser  beam 
emanating  from  a  circular  or  annular  aperture  in  the  presence  of  turbulence, 

(ii)  to  estimate  the  improvement  in  mean  focal  intensity  obtainable  by  adaptive 
control  of  segmented  apertures. 

The  programs  are  based  on  expressions  taken  or  developed  from  Lutomirksi* .  It 
must  be  emphasised  that  the  theory  is  far  from  exhaustive  and  in  particular  is 
only  applicable  to  low-power  laser  sources.  Hence,  as  is  often  the  case,  the 
accuracy  of  the  numerical  solution  exceeds  that  merited  by  the  theory. 

The  programs  provide  close  agreement  with  Lutomirski's  graphs  and  with 
exact  expressions  available  for  free  space  propagation.  The  first  program  loses 
accuracy  as  the  longitudinal  distance  (characterised  by  8  )  from  the  focus 
increases  but  this  could  be  overcome  by  implementing  a  more  appropriate  integra¬ 
tion  method  (Filon's  would  be  adequate)  for  the  inner  integral  in  equation  (II). 
The  second  program  only  predicts  the  mean  intensity  at  the  focal  point. 
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Appendix  A 

USING  PROGRAM  M2 1 LASER 

This  program  evaluates  the  integral  in  (12)  using  QAOZA  of  the  Harwell 
Subroutine  Library.  Jq  is  obtained  via  the  ICL  subroutine  F4J0  while  the  cal¬ 
culation  of  Ms(x,z)  is  described  in  section  3.2  and  M(6,B,x)  is  covered  in 
section  3.1  and  extended  to  include  annular  apertures  (£  #0)  in  section  4. 

The  data  consists  of  (NRUN  +4)  cards. 

(1)  NRUN:  gives  the  number  of  different  values  of  a  which  will  be  specified. 
NRUN  *  0  stops  the  program. 

(2)  z/zc  ,  D/Lq  ,  ELIM: 

z/z£  =  0  for  free  space.  ELIM  determines  whether  the  x-range  will  be 

split  into  two  (0  ^  x  ^  0.25  ,  0.25  ^  x  ^  1)  .  ELIM  *  0  forces  this  to 

occur,  otherwise  ELIM  =4  is  recommended  so  that  the  range  will  be  split 

if  M  is  decreasing  quickly  (see  section  3.3). 
s 

(3)  6,  5,  £ 

(4)  ERR:  The  absolute  accuracy  required  in  the  value  of  the  integral * 

ERR  =  0.0001  is  recommended  but  it  can  of  course  be  reduced  if  the  inte¬ 
gral  unexpectedly  turns  out  to  be  small. 

(5)  a  :  The  program  will  calculate  the  integral  and  then  ask  for  the  next  value 
of  a  until  NRUN  values  have  been  processed  after  which  a  new  value  for 
NRUN  is  expected. 

The  output  consists  of  a  single  line 

VALUE  OF  INTEGRAL  =  X.XXXXEXXX  ESTIMATED  ERROR  -  X.XXXXEXXX  XXX  INTEGRAND 
EVALUATIONS  . 


The  Strehl  ratio  can  then  be  obtained  by  dividing  by  expression  (27)  if  a  *  0 

or  by  (31)  if  a  ¥=  0  but  6  =  0  *  £  -  0  .  Otherwise,  the  intensity  in  vacuo 

has  to  be  evaluated  numerically  by  running  the  program  with  z/z  ■  0  . 

c 

All  integrals  are  evaluated  by  Simpson’s  method,  the  maximum  interval 
width  (m.i.w.)  having  been  chosen  empirically.  M(6,0,£,x)  is  found  from  equa¬ 
tion  (11),  the  range  of  the  inner  integral  being  modified  for  annular  apertures 
as  described  in  section  4.  The  m.i.w.  are  specified  in  DATA  statements  and  are 
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Appendix  A 


Likewise,  the  m.i.w.  for  the  integrals  in  (25)  involved  in  the  determination  of 
Ms(x,z)  when  Dx/Lq  >  3  ,  are  specified  in  DATA  statements. 

HODU  =  -  =10  ;  N  =  50  (hence  Aw  *  0.02wj) 

max 

Finally,  QAOZA  needs  to  know,  when  evaluating  (12),  the  smallest  permis~ 

sible  interval.  Ax  .  and  this  has  been  set  at  DXMIN  *  0.007  .  QAOZA  begins 

min 

by  dividing  the  range  into  16  equal  intervals  and  then  successively  halves  the 
step  locally  to  obtain  the  desired  accuracy.  Thus  when  the  interval  (0,1)  is 
treated,  a  maximum  of  127  integrand  evaluations  will  be  performed 
(Ax  =  0.0078125)  and  when  it  is  divided  into  (0,0.25)  and  (0.25,1)  the 
largest  number  possible  is  32  +  63  (Ax  *  0.01171875)  =  95  ,  where  we  have  used 
the  fact  that  the  integrand  is  zero  at  x  *  0  and  x  =  1  . 

When  NRUN  >  1  we  need  to  recalculate  only  Jq  in  equation  (12)  when 
a  changes.  We  therefore  store  M((S,6,£,x)  and  M  (x,z)  (in  arrays  AXMDB, 

AEMS)  and  record  which  values  of  x  have  occurred  (in  array  LST)  in  order  to 
avoid  unnecessary  work.  Thus  the  time  taken  to  treat  NRUN  values  of  a  is 
much  less  than  NRUN  times  the  time  for  one  a  .  However,  should  one  wish  to 
reduce  DXMIN  then  the  size  of  these  arrays  must  be  increased  accordingly. 
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Appendix  B 

USING  PROGRAM  M21QPTSR 

This  program  evaluates  the  integral  in  equation  (34),  M  being  calcula- 

s 

ted  as  described  in  section  3.2  (with  HODU  =  10  and  N  =  50  when  applicable, 
as  in  Appendix  A).  The  data  consists  of  three  records  as  follows 

(1)  f/zc  »  :  ^zc  =  ®  f°r  ^ree  space  and  f / Zq  <  0  stops  the  program. 

(2)  HMAX,  HF,  DTH:  recommended  values  0.1,  1.0,  0.05 

(3)  6,  NS:  NS  is  the  number  of  equal  sectors  comprising  the  annulus 

(ze  0  *  2ir/NS)  . 


Each  range  of  integration  is  divided  into  an  even  number  of  intervals  and 
Simpson's  rule  applied.  The  maximum  interval  widths  are 


A0 


max 


DTH 


Ax 


max 


mm 


HMAX 


The  output  from  the  program  consists  of  a  single  line 

STHEHL  RATIO  FOR  SECTOR  *  X.XXXXEXXX 

where  the  answer  is  obtained  by  dividing  (34)  by  (35). 
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Appendix  C 

THE  PROBABILITY  DENSITY  FUNCTION  FOR  THE  SEPARATION 
OF  TWO  POINTS  PLACED  RANDOMLY  WITHIN  AN  ANNULUS 

Let  the  outer  diameter  be  D  and  the  inner  diameter  be  £D.  The  p.d.f*  p(Dx) 
is  related  to  the  area  of  overlap,  A  ,  of  two  identical  annuli  whose  centres  are 
separated  by  Dx  . 


A  =  Aj  +  A^  ~  2A2 


(C-l) 


A,  is  the  area  of  overlap  of  the  outer  circles 


Aj  =  |D2(cos  5x  -  x(l  -  x2)^) 


(C-2) 


A 2  is  the  area  enclosed  between  an  inner  circle  and  the  outer  circle  which 
is  not  concentric  with  it. 


2A„  =  |D2  £2 


jcos  *b  -  b(  1  -  b2)^j  +  £cos  *a  -  a(l  ~  a2)^J 


(C-3) 


where 


a  =  l^i-+x 
4x 


(i  -  r 


(C-4) 


A^  is  the  area  of  overlap  of  the  inner  circles 


A3  -  I (5D)2 jcos-1 (x/5)  -  (x/5)^l  -  (x/C)2 
The  'normalised'  p.d.f.  is  then 


(C-5) 


P(Dx)  •  ^yA 
ttD 


(C-6) 


so  that 


P(Dx)dx  -  (1  -  £2)2 


(C-7) 
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Note  that  A^  =  0  when  x  >  5  while  from  (C-4)  one  can  deduce  that 


a,  [b| 

> 1  when 

X 

>  id 

+  O  or  X  <  1(1  -  o 

giving 

^2 

■  0  when 

X 

>  10 

+  O  ; 

^2 

-  K2d2 

when 

X 

<  id  -  O 

The  latter  implies  that  the  inner  circle  is  wholly  contained  within  the  other 
outer  circle  when  x  <  J(1  -  O  (as  in  Fig  3c). 
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Table  1 


SHOWING  THE  IMPROVEMENT  IN  STREHL  RATIO  PREDICTED 


BY  OPTIMISING  THE  SEGMENTED  APERTURE  OF  FIG  4b 


c2 

n 

f/zc 

n2 

^opt 

n annul us 

io-12 

3000 

0.0882 

0.0864 

0.0876 

0.0059 

uf13 

300 

0.588 

0.580 

0.585 

0.0874 

io'14 

30 

0.941 

0.940 

0.941 

0.600 

io-15 

3 

0.994 

0.994 

0.994 

0.945 

.o-16 

0.3 

0.9994 

0.9994 

0.9994 

0.995 

2  — 

C  is  in  units  of  m  J  .  Columns  3  and  4  were  obtained  with  M210PTSR,  the 
n 

fifth  column  being  derived  from  them  using  (37).  M21LASER  produced  column  6. 
Hence  in  this  example  the  optimally  adjusted  segmented  aperture  performs  as 
well  as  the  full  aperture  would  in  ten  times  weaker  turbulence. 


REFERENCES 


Author 

R.F.  Lutomirski 

R.F.  Lutomirski 
H.T.  Yura 

Y.L.  Luke 

D.M.  Gilbey 


Title,  etc 

Propagation  of  focused  truncated  laser  beams  in  the 
atmosphere. 

In  Optical  propagation  in  the  atmosphere.  AGARD  CP- 183 
Paper  22  (1976) 

Wave  structure  function  and  mutual  coherence  function  of 
an  optical  wave  in  a  turbulent  atmosphere. 

J.  Opt.  Soc.  Amer.y  6J_,  482-491  (1971) 

Integrals  of  Bessel  functions. 

In  Handbook  of  mathematical  functions .  ed.  M.  Abramowitz 
and  I. A.  Stegun,  Ch  1 1 . 

Dover,  New  York  (1972) 

Unpublished  RAE  work 

AF~  rot  f.  CzySARILY 
;v  PUBLIC 
.-!AL  CBGAM'i.'iTIONS 


y.  .  :  \Z.  . 

Of,  T.o  C0M.\ 


M-C  7903 


Figs  1&2 


Fig  1  Basic  geometry  showing  two  typical  points  in  the  aperture 
and  the  point,  P  at  which  the  intensity  is  required 


Fig  2  Two  alternative  parametrisations  of  the  domain  of  the 
inner  integral  in  equation  (4) 


enclosed  between  two  identical  overlapping  annuli 


Fig  4 


Fig  4  Segmented  apertures 


t  (petres) 


Fig  7 


(*)  No  tor balance  ;  3  •  0.5 

2  -15  -2/3 

(b)  C#  •  3  x  10  cm  ;  J  ■  0.5 


Comparison  of  on-axis  Intensities  at  x  =  10. 6v  for 
D  «  30  cm  focussed  at  0.5  km  for  a  uniform  aperture 
distribution  (af  Fig  L8a). 
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